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<^ Quantum simulation is an important application of future quantum computers 

with applications in quantum chemistry, condensed matter, and beyond. Quantum 
simulation of fermionic systems presents a specific challenge. The Jordan- Wigner 
transformation allows for representation of a fermionic operator by 0(n) qubit op- 
erations. Here we develop an alternative method of simulating fermions with qubits, 

a 

first proposed by Bravyi and Kitaev [S. B. Bravyi, A.Yu. Kitaev, Annals of Physics 

O" 1 

298, 210-226 (2002)], that reduces the simulation cost to O(logn) qubit operations 
for one fermionic operation. We apply this new Bravyi-Kitaev transformation to the 
task of simulating quantum chemical Hamiltonians, and give a detailed example for 

On 

I/") the simplest possible case of molecular hydrogen in a minimal basis. We show that 

00 

O the quantum circuit for simulating a single Trotter time-step of the Bravyi-Kitaev de- 

<N 

y—i rived Hamiltonian for H2 requires fewer gate applications than the equivalent circuit 

derived from the Jordan- Wigner transformation. Since the scaling of the Bravyi- 

X 

Kitaev method is asymptotically better than the Jordan- Wigner method, this result 
for molecular hydrogen in a minimal basis demonstrates the superior efficiency of 
the Bravyi-Kitaev method for all quantum computations of electronic structure. 



I. INTRODUCTION 

In his seminal article that anticipated the field of quantum information, Feynman argued 
that simulating quantum systems on classical computers takes an amount of time that scales 
exponentially with the size of the system, while the cost of quantum simulations can scale in 
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polynomial time with system size pp. This possibility may offer a path forward for computa- 
tional chemistry [2113] . A quantum simulation algorithm for quantum chemical Hamiltonians 
enables the efficient calculation of properties such as energy spectra [3] , reaction rates (3J [5] , 
correlation functions [6], and molecular properties [7] for molecules larger than those that 
are currently accessible through classical calculations. 

Quantum simulation of electronic structure requires a representation of fermions by sys- 
tems of qubits. Significant progress has been made on efficient quantum simulation of 
fermions. In 1997, Abrams and Lloyd proposed a simulation scheme for fermions hopping 
on a lattice [8]. In 2002, Somma et al. used the Jordan- Wigner to generalize the simulation 
scheme proposed by Abrams and Lloyd [SI HO]- The Jordan- Wigner transformation has since 
been used to outline a scalable quantum algorithm for the simulation of molecular electron 
dynamics, and to design an explicit quantum circuit for simulating a Trotter time-step of 
the molecular electronic Hamiltonian for H 2 in a minimal basis [31 [11] . Further refinements 
of the Jordan- Wigner construction were made by Verstrate and Cirac [12] and by Bravyi 
and Kitaev [13J. From the point of view of fundamental physics, such constructions can be 
regarded as giving a negative answer to the question of whether fundamental fermi fields 
are required to explain observed fermionic degrees of freedom [14]. Practically speaking, 
such constructions show that quantum computation of electronic structure does not suffer 
from an analog of the sign problem; that is, fermion antisymmetry represents no significant 
obstacle to efficient algorithms. 

Theoretical progress in quantum simulation has been accompanied by experimental suc- 
cesses. In 2010, Lanyon et al. calculated the energy spectrum of a hydrogen molecule using 
an optical quantum computer [15] . For a review of photonic quantum simulators, see [T6] . 
Du et al. repeated this result to higher precision with NMR shortly thereafter [17] . Digital 
quantum simulations of the kind considered in the present paper have been implemented 
in ion traps using up to 100 gates and 6 qubits [IB]. The progress of trapped ion quantum 
simulation is detailed in [T9] . 

Quantum computation of electronic structure has been the subject of simulation stud- 
ies [HIEO] and has been extended to cover relativistic systems [21]. The history of calculations 
in quantum chemistry provides a useful sequence of problems reaching from calculations that 
can be performed on experimental quantum computers today to calculations at the present 
research frontier [22] • Despite these promising results, the scaling of the number of gates 
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required by the algorithm outlined in [3l [TT] remains challenging. It is a subject of active 
research to find improvements to the (polynomial) scaling of the cost of the algorithm de- 
scribed in 0, [EE] • Several improvements are described in [23] , and the techniques of that 
work could be combined with those of the present paper to further reduce the resource 
requirements. 




FIG. 1: A simulation scheme first encodes fermionic states in qubits, then acts with the qubit 
operator representing the fermionic operator (obtained by the associated transformation), then 
inverts the encoding to obtain the resultant fermionic state. The criterion for a successful simulation 
scheme is that this procedure reproduces the action of the fermionic operator, i.e. that Path 1 is 
equivalent to Path 2, for all basis states — in other words, that this diagram commutes. 

A fermionic simulation scheme can be broken into two pieces: first, to map occupation 
number basis vectors to states of qubits; and second, to represent the fermionic creation 
and annihilation operators in terms of operations on qubits in a way that preserves the 
fermionic anti-commutation relations, as illustrated in Figure [TJ Previous simulation algo- 
rithms have used a straightforward mapping of fermionic occupation number basis states to 
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qubit states that was originally defined by Zanardi in the context of entanglement [31 [9j |2l] . 
The Jordan- Wigner transformation is then used to write the electronic Hamiltonian as a 
sum over products of Pauli spin operators acting on the qubits of the quantum computer. 
Subsequently the Hamiltonian terms hk, where H = J2k^k, are converted into the unitary 
gates that are the corresponding time evolution operators. Even though the hk do not 
necessarily commute, their sequential execution on a quantum computer can be made to ap- 
proximate the unitary propagator e~ tHt through a Trotter decomposition [25H28] . Finally, 
the iterative phase estimation algorithm (IPEA) is used to approximate the eigenvalue of 
an input eigenstate [31 EH [28] ■ 

In this paper we treat the Trotterization process and IPEA as standard procedures. We 
develop the Bravyi-Kitaev basis and Bravyi-Kitaev transformation, both named after the 
authors who first proposed such a scheme [13], which provide a more efficient mapping 
between electronic Hamiltonians and qubit Hamiltonians. While the occupation number 
basis and the Jordan- Wigner transformation allow for the representation of a single fermionic 
creation or annihilation operator by 0(n) qubit operations, the Bravyi-Kitaev basis and 
transformation require only O(logra) qubit operations to represent one fermionic operator. It 
is worth noting that Bravyi and Kitaev were concerned with exploring the power of fermions 
as the basic hardware units of a quantum computer, rather than with the simulation of 
fermions by qubits [13] . However, understanding how the structure of fermionic systems can 
be employed to process information helps us understand how standard quantum information 
procedures can be used to simulate the structure of fermionic systems. We work out a 
detailed application of the Bravyi-Kitaev transformation to the operators that appear in 
quantum chemical Hamiltonians, providing a new way of mapping electronic Hamiltonians 
to qubit Hamiltonians. We also give explicit Pauli decompositions of the qubit operators 
derived from this new transformation for the quantum chemical Hamiltonian for H 2 in a 
minimal basis. We show that the quantum circuit for simulating a single first-order Trotter 
time-step of the Bravyi-Kitaev minimal basis molecular hydrogen Hamiltonian requires 30 
single-qubit gates and 44 CNOT gates, as compared to 46 single-qubit gates and 36 CNOT 
gates for the Jordan- Wigner Hamiltonian derived in [11] . Finally, we show that a chemical- 
precision estimate of the ground state eigenvalue of the Bravyi-Kitaev Hamiltonian can be 
obtained in 3 first-order Trotter steps, with a total cost of 222 gates, while the Jordan- 
Wigner Hamiltonian requires 4 first-order Trotter steps for a total of 328 gates. Since the 
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Bravyi-Kitaev transformation is known to be asymptotically more efficient, this result for the 
simplest possible case of molecular hydrogen in a minimal basis demonstrates the superior 
efficiency of the Bravyi-Kitaev method for all molecular quantum simulations. 

In Section [TT] we will review basic quantum chemistry in second quantized form as well as 
the Jordan Wigner transformation. In Section |III| we discuss alternatives to the occupation 
number basis, including the Bravyi-Kitaev basis, which we go on to describe in detail in 



Section |IV| In Section [V] we present the Bravyi-Kitaev transformation, which allows us to 

we 



represent creation and annihilation operators in the Bravyi-Kitaev basis. In Section VI 



compute the products of these operators that occur in electronic structure Hamiltonians. In 



Section VII we compute the molecular electronic structure Hamiltonian of H2 in a minimal 



basis using the Bravyi-Kitaev basis and transformation. In Section VIII we make an explicit 
comparison between the Bravyi-Kitaev transformation and the Jordan Wigner transforma- 
tion by simulating the Trotterization procedure. We close the paper with some conclusions 
about the utility of the Bravyi-Kitaev transformation. 



II. BACKGROUND 



A. Fermionic systems and second quantization 

We may describe fermionic systems using the formalism of second quantization, in which 
n single-particle states can be either empty or occupied by a spinless fermionic particle. In 
the context of quantum chemistry these n states represent spin orbitals, ideally one-electron 
energy eigenfunctions and often molecular orbitals found by the Hartree-Fock method [29l 
[30] . We consider a subspace of the full Fock space which is spanned by 2 n electronic basis 
states \ f n -i ... fo), where fj G {0,1} is the occupation number of orbital j (restricted 
to these values due to the Pauli exclusion principle). This is called the occupation number 
basis. 

Any interaction of a fermionic system can be expressed in terms of products of the 
creation and annihilation operators aj and aj, for j G {0, . . . , n— 1}. Due to the exchange 
anti-symmetry of fermions, the action of aj or aj introduces a phase to the electronic basis 
state that depends on the occupancy of all orbitals with index less than j in the occupation 
number representation. (One can choose instead to define these operators so that it is the 
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occupation of orbitals with index greater than j that determines the phase — the ordering 
of orbitals is arbitrary.) These operators act on occupation number basis vectors as follows: 

-i) E -' /s l/„-i ••• I j , , i ./• , ••• /o); (i) 

= 0; (2) 
;_!)££/. 1/^ ... ./•,,()/, , ... / ); (3) 

= 0. (4) 

The canonical fermionic anti-commutation relations enforce the exchange anti-symmetry: 

[aj, a k ] + = 0, [a], a{} + = 0, [aj, a\} + = 6 jk l, (5) 

where the anti-commutator of operators A and B is defined by [A, B] + = Ai? + 

The molecular electronic Hamiltonian of interest in the electronic structure problem is: 

H = ^ hjj aja^ + - ^ /^ijjfci a|ata fc a z . (6) 

The coefficients hy and /iyj.; are one- and two-electron overlap integrals, which can be 
precomputed classically and input to the quantum simulation as parameters [21 [HJ |2H] ■ 



As an application of the techniques presented in this paper (Section VII), we treat molec- 
ular hydrogen in a minimal basis. Thus, we construct two spatial molecular orbitals by 
taking linear combinations of the localized atomic spatial wavefunctions: ip g = ipni + ipH2 
and ip u = ipHi — ipH2- Here the subscripts g and u stand for the German words gerade and 
ungerade — even and odd. In general one must take a Slater determinant to determine 
the correctly anti-symmetrized wavefunctions of the fermionic system, but in this case we 
can guess them by inspection. The form of the spatial wavefunctions is determined by the 
choice of basis set. STO-3G is a commonly used Gaussian basis set — for further details 
see [291130]. 

Molecular spin orbitals are formed by taking the product of these two molecular spatial 
orbitals with one of two orthogonal spin functions, \a) and |/3). Thus, the four molecular 
spin orbitals in our model of the hydrogen molecule (which correspond to the operators a 



(th 



arc: 



IXo) = |V> fl >|a), l*i) = hM0>, |X2> = hk.)|a>, IXs) = (7) 
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In the next section we will review the occupation number basis and the Jordan-Wigner 
transformation, which together have been established as a standard method for mapping 
fermionic systems to quantum computers [31 El EH EE] • 

B. The Jordan-Wigner transformation 

The form of electronic occupation number basis vectors suggests the following identifica- 
tion between electronic basis states on the left and states of our quantum computer [21]: 

|/n-i • • • h fo) -»> k-i) • • • ® |<?i> ® |?o), /j = e {0, 1}. (8) 

That is, we let the state of each qubit \qj) store fj, the occupation number of orbital j. We 
refer to this method of encoding fermionic states as the occupation number basis for qubits. 
The next step is to map fermionic creation and annihilation operators onto operators on 
qubits. 

We can form one-qubit creation and annihilation operators, Q + and Q~, that act on 
qubits of our quantum computer as follows: 

g+|o) = |i), g+|i) = o, <r|i> = |o>, q-|o) = o. (9) 

We could proceed by following the standard recipe for turning p-qubit quantum gates into 
operators acting on an n-qubit quantum computer (n > p) by taking the tensor product of 
the gates acting on the target qubits with the identity acting on the other (n — p) qubits. 
However, it is easy to show that the qubit creation and annihilation operators formed in 
this way do not obey the fermionic anti-commutation relations. 

Expressing the qubit creation and annihilation operators in terms of Pauli matrices sug- 
gests a way forward: 

Q+ = |1)(0| = \{a x - io"), Q- = |0)(1| = ~K + io*). (10) 

The mutual anti-commutation of the three Pauli matrices allows us to recognize that Q ± 
anti-commutes with a z . Thus if we represent the action of al or aj by acting with and 
with o z on all qubits with index less than j, our qubit operators will obey the fermionic anti- 
commutation relations. Put differently, the states of our quantum computer will acquire the 
same phases under the action of our qubit operator as do the electronic basis states under 
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the action of the corresponding creation or annihilation operator. The effect of the string 
of o z gates is to introduce the required phase change of —1 if the parity of the set of qubits 
with index less than j is 1 (odd), and to do nothing if the parity is (even), where the parity 
of a set of qubits is just the sum (mod 2) of the numbers that represent the states they are 
in. 

We can then completely represent the fermionic creation and annihilation operators in 
terms of basic qubit gates as follows: 

at = I®"-''" 1 <g> Q+ <g> [o* 9j ], aj = l® n - j - 1 ®Q~ ®[a z ® j }. (11) 

A more compact notation, of which we will make extensive use throughout this paper, is: 

a] = Qf <g> Zf_ x = \{Xj <8> Zf_ x - iYj <g> Zjt x )\ (12) 

aj ee QJ <g> Zjli = \{Xj ® ZJL X + iYj <g> Z^), (13) 

where: 

Z? ee at ®al 1 ®---a{® a*, (14) 

and where it is assumed that any qubit not explicitly operated on is acted on by the identity. 
The operator Z~* is a "parity operator" with eigenvalues ±1, corresponding to eigenstates 
for which the subset of bits with index less than or equal to % has even or odd parity, 
respectively. 

The above correspondence, a mapping of interacting fermions to spins, is the Jordan- 
Wigner transformation [31 [101 IH El] . Jordan and Wigner introduced this transformation 
in 1928 in the context of ID lattice models, but it has since been applied to quantum 
simulation of fermions [31 19TTTT] . The problem with this method is that as a consequence of 
the non-locality of the parity operator Z^, the number of extra qubit operations required 
to simulate a single fermionic operator scales as 0(n). In the next section we consider two 
alternatives to the occupation number basis that were suggested by Bravyi and Kitaev [13]. 

III. ALTERNATIVES TO THE OCCUPATION NUMBER BASIS 

A. The parity basis 

The extra qubit operations required to simulate one fermionic operator when using the 
Jordan- Wigner method result from operating with a z on all qubits with index less than j. 
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This task could be accomplished by a single application of a z if instead of using qubit j 
to store fj, we used qubit j to store the parity of all occupied orbitals up to orbital j [T3] . 
That is, we could let qubit j store pj = J2l=o fs- (Throughout this paper, all sums of binary 
variables are taken mod 2). We follow [13] and call this encoding of fermionic states in 
qubit states the parity basis. 

It is useful to define the transformations between bases we will consider in terms of maps 
between bit strings. For all the transformations we consider, which involve only sums of bits 
mod 2, it is possible to represent their action by matrices acting on the vector of bit values 
corresponding to a given logical basis state. For example, the occupation number basis state 
I/7 •• • /1/0) is equivalent to the following vector: 



C/Vj • • • > .fi> fo) T 

In terms of these vectors the map to the parity basis is given by: 

Pi = '^^[ K n}ij fji 



(15) 



(16) 



where n is the number of orbitals. 7r n is the (n x n) matrix defined below. Note that we 
index the matrix tt„ from the lower right corner, for consistency with our orbital numbering 
scheme. 



7T 



so that TT r , 



(17) 



1 i < j 
i>j 

- • • I 

For example, to change the occupation number basis state 1 10100111) into its corresponding 
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parity basis state 1 10011101), we act with the matrix 7r§ on the appropriate bit string: 
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(18) 



With this understanding of the parity basis transformation, we can now derive the transfor- 
mation that maps fermionic operators into operators in the parity basis. Since the parity of 
the set of orbitals with index less than j is what determines whether the action of a*- intro- 
duces a phase of —1, operating with a z on qubit (j — 1) alone will introduce the necessary 
phase to the corresponding qubit state in the parity basis. 

However, unlike the Jordan- Wigner transformation, we cannot represent the creation or 
annihilation of a particle in orbital j by simply operating with Q ± on qubit j, because in the 
parity basis qubit j does not store the occupation of orbital j, but the parity of all orbitals 
with index less than or equal to j. Thus whether we need to act with Q + or Q~ on qubit 
j depends on qubit (j — 1). If qubit (j — 1) is in the state |0), then qubit j will accurately 
reflect the occupation of orbital j, and simulating aj will require acting on qubit j with Q + , 
as before. But if qubit (j — 1) is in the state |1), then qubit j will have inverted parity 
compared to the occupation of orbital j, and we will instead need to act with Q~ on qubit 
j to simulate a] (and vice versa for the annihilation operator). 

The operator equivalent to Q ± in the parity basis is therefore a two-qubit operator acting 
on qubits j and j — 1: 

Vf = Qf <g> |0> <0|^- ! - QJ g) llXl^x = 1 -{X 3 ® T iTj). (19) 

Additionally, creating or annihilating a particle in orbital j changes the parity data that must 
be stored by all qubits with index greater than j. Thus we must update the cumulative sums 
Pk for k > j by applying a x to all qubits \pk), k > j [13]. The representations of the creation 
and annihilation operators in the parity basis are then: 
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a) = Xj- +1 <g> P/ = f (Jfl^ ® ® - zXj^ (8) Fj); 



(20) 
(21) 



where: 



(22) 



This is the equivalent of the Jordan- Wigner transformation for the parity basis. The operator 
Xf~ is the "update operator", which updates all qubits that store a partial sum including 
orbital (i — 1) when the occupation number of that orbital changes. It is straightforward 
to verify that these mappings satisfy the fermionic anti-commutation relations. But to 
simulate fermionic operators in the parity basis, we have traded the trailing string of o z 
gates required by the Jordan- Wigner transformation for a leading string of a x gates whose 
length also scales as O (n) , and we have not improved on the efficiency of the Jordan- Wigner 
simulation procedure. In the next section, we explore a third possibility. 



Two kinds of information are required to simulate fermionic operators with qubits: the 
occupation of the target orbital, and the parity of the set of orbitals with index less than 
the target orbital. The previous two approaches are dual in the way that they store this 
information. With the occupation number basis and its associated Jordan- Wigner transfor- 
mation, the occupation information is stored locally but the parity information is non-local, 
whereas in the parity basis method and its corresponding operator transformation, the parity 
information is stored locally but the occupation information is non-local. 

The Bravyi-Kitaev basis is a middle ground. That is, it balances the locality of occupation 
and parity information for improved simulation efficiency. The general form of such a scheme 
must be to use qubits \bj) to store partial sums J2 l s =k fs °f occupation numbers according 
to some algorithm. For ease of explanation, in the exposition that follows, when we write 
that a qubit "stores a set of orbitals" , what is meant is that the qubit stores the parity of 
the set of occupation numbers corresponding to that set of orbitals. 

Bravyi and Kitaev's encoding has an elegant binary grouping structure [13]. In this 
scheme, qubits store the parity of a set of 2 X orbitals, where x > 0. A qubit of index j 



B. The Bravyi-Kitaev basis 
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always stores orbital j. For even values of j, this is the only orbital that it stores, but for 
odd values of j, it also stores a certain set of adjacent orbitals with index less than j. Just 
as with the parity basis transformation, this encoding can be symbolized in a matrix j3 n that 
acts on bit string vectors corresponding to occupation number basis vectors of length n to 
transform them to the corresponding Bravyi-Kitaev-encoded bit strings (again, all additions 
done mod 2). In terms of these vectors, the map from the occupation number basis to the 
Bravyi-Kitaev basis is: 

h = EiPnk fj, (23) 
3 

where the matrix (3 n is given in Figure [2] below. 



1 














P2* 
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P; 



FIG. 2: The matrix /3 n that transforms occupation number basis vectors of length n into the 
Bravyi-Kitaev basis. /3i is a (1 x 1) matrix with a single entry of 1. Subsequent iterations of the 
matrix that act on occupation number basis vectors of length 2 X are constructed by taking ligi/^-i 
and then filling in the top row of the first quadrant of this matrix with l's. j3 n for 2 X < n < 2 X+1 
is just the (n x n) segment of f3 2 x+i that includes 60 through b n —\. The recursion pattern for the 
inverse transformation matrix is also shown. An entry of 1 in row bi, column fj means that bi is a 
partial sum including fj. 

For example, to change the occupation number basis state 1 10100111) into its correspond- 
ing Bravyi-Kitaev basis state 1 10101101), we act with the matrix /3g on the appropriate bit 
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string vector: 
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(24) 



This encoding strikes a balance between the occupation number basis and the parity 
basis methods. The parity of occupied orbitals up to orbital j is no longer stored in a single 
qubit, but the Bravyi-Kitaev encoding stores the parity of orbitals with index less than j in 
a few partial sums whose number scales as O(logj) < O(logn) [15] . Likewise, we no longer 
need to update all the qubits with index greater than j, but only those that store partial 
sums which include occupation number j. Each occupation number enters an additional 
partial sum only if the number of single particle states n is doubled, and so the overall cost 
of simulating a single fermionic operator with qubits scales as O(logn) [13J. 

Given this encoding, we need to determine — for an arbitrary index j — which qubits in 
the Bravyi-Kitaev basis store the parity of all orbitals with index less than j, which qubits 
store a partial sum including orbital j, and which qubits determine whether qubit j has the 
same parity or inverted parity with respect to orbital j. These sets of indices will allow us to 
explicitly construct the fermionic creation and annihilation operators in the Bravyi-Kitaev 
basis. In the next section, we define these sets of qubit indices. 

IV. SETS OF QUBITS RELEVANT TO THE BRAVYI-KITAEV BASIS 

In this section we define the sets of qubits that are involved in the Bravyi-Kitaev trans- 
formation. These are the parity set (the qubits in the Bravyi-Kitaev basis that store the 
parity of all orbitals with index less than j), the update set (the qubits that store a partial 
sum including orbital j), and the flip set (the qubits that determine whether qubit j has the 
same parity as orbital j). 
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A. The parity set 

For an arbitrary index j, we would like to know which set of qubits in the Bravyi-Kitaev 
basis tells us whether or not the state of the quantum computer needs to acquire a phase 
change of —1 under the action of a creation or annihilation operator acting on orbital j. 
The parity of this set of qubits has the same parity as the set of orbitals with index less 
than j, and so we will call this set of qubit indices the "parity set" of index j, or P(j). 
To determine the elements of P(j), we consider the transformation from the Bravyi-Kitaev 



basis to the parity basis. From equation (16) we know that Pi = J2j[^n]ij fj- Given the 
inverse transformation matrix f3~ l , it is also true that: 



and hence: 



Si = Et/^U h, 

k 



Pi = Ej[^n}ij (Efc[/VU h) 



(25) 



(26) 
(27) 



The matrix n n ^n l ls the transformation matrix from the Bravyi-Kitaev basis to the parity 
basis. Therefore, the nonzero entries to the right of the main diagonal in row i of the matrix 
7r n /3~^ give the indices of qubits in the Bravyi-Kitaev basis that can be used to compute the 
cumulative parity of orbitals with index less than i. An entry of 1 in row i, column j of 
7r n /3~^ (where j < i, i.e. to the right of the main diagonal by our numbering) indicates that 
j e P(i): 
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B. The update set 

For arbitrary j, we define the set of qubits (other than qubit j) that must be updated when 
the occupation of orbital j changes. We call this set the "update set" of index j, or U(j). 
This is the set of qubits in the Bravyi-Kitaev basis that store a partial sum including orbital 
j. Any Bravyi-Kitaev qubit that stores a partial sum that includes occupation number j is 
in U(j). Since even indexed qubits store only the occupation of the corresponding orbital, 
update sets contain only odd indices. It is straightforward to determine the elements of U (?) 
from the transformation matrix (3 n that maps bit strings in the occupation number basis 
to the Bravyi-Kitaev basis. The columns of this transformation matrix show which qubits 
in the Bravyi-Kitaev basis store a particular orbital, and so the nonzero entries in column 
j above the main diagonal determine the qubits other than qubit j that must be updated 
when the occupancy of orbital j changes. These are the elements of the update set. 
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(29) 



It should be clear that update sets depend on the size of the basis used. For example, if 
16 basis functions were used instead of the 8 used in the example above, all the update sets 
other than U(7) would also include index 15. 

C. The flip set 



For arbitrary j, we need to know what set of Bravyi-Kitaev qubits determines whether 
qubit j has the same parity or inverted parity with respect to orbital j. We will call this 
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set of Bravyi-Kitaev qubits the "flip set" of j, or F(j), because this set is responsible for 
whether bj has flipped parity with respect to /_,-. This is the set that stores the parity of 
occupation numbers other than fj in the sum bj. Since even-indexed qubits store only the 
orbital with the same index, the flip set of even indices is always the empty set. One can 
determine the elements of F(j) by looking at the inverse transformation matrix that 
maps bit strings in the Bravyi-Kitaev basis to the occupation number basis. The columns 
with nonzero entries to the right of the main diagonal in row % of this inverse transformation 
matrix give the indices of the Bravyi-Kitaev qubits that together store the same set of 
orbitals as is stored by These are the elements of the flip set. 
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(30) 



With these sets defined, we can derive the mapping from fermionic operators to qubit op- 
erators that is the equivalent of the Jordan- Wigner transformation in the Bravyi-Kitaev 
basis. 



V. THE BRAVYI-KITAEV TRANSFORMATION 



In this section we will give an explicit prescription, in terms of Pauli matrices, for repre- 
senting the creation and annihilation operators that act on the Bravyi-Kitaev basis states. 
Operating in this basis requires that we find the analogues to the qubit creation and anni- 
hilation operators {Q^ in the occupation number basis, V ± in the parity basis) as well as 
the parity operator, Z^, and the update operator, Xf~, in the Bravyi-Kitaev basis. We will 
first define some notation. 
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For our purposes it is the parity of subsets of orbitals or qubits that matters, not the 
individual occupation numbers or states of the qubits in the set. Thus, it is useful to define 
operators that project onto the subspace of the Hilbert space of the entire computer for 
which the subset of qubits with indices in S has the parity selected for by the operator (even 
for Es, odd for O5). We can express these operators in terms of Pauli matrices as follows: 

E s =^{l + Z s ), O s = ^(l-Z s ), (31) 

where Z$ is shorthand for the a z gate applied to all qubits in S. With this notation 
established, we will next write equations for the qubit operators in the Bravyi-Kitaev basis 
that represent creation and annihilation operators acting on orbital j. To begin we will 
consider the case for which j is even, because this will allow us to build intuition for the 
more difficult case for which j is odd. 

A. Representing in the Bravyi-Kitaev basis for j even 

In the case that j is even, we should act with Q on qubit j, just as for the Jordan-Wigner 
transformation, because the Bravyi-Kitaev encoding stores orbitals with j = (mod 2) in 
the qubit with the same index. There are then two additional tasks that dictate how to 
represent the fermionic operators in the Bravyi-Kitaev basis: determining the parity of 
occupied orbitals with index less than j, and updating qubits with index greater than j that 
store a partial sum that includes occupation number j. 

The parity of the set of qubits in P(j) is equal to that of the set of orbitals with index less 
than j. By analogy with the Jordan-Wigner transformation, we act with a z on all qubits 
with indices in P(j), that is, we apply the operator Zp(jy The number of qubits in P(j) 
scales as O(logj) < O(logn) [T3] . 

Secondly, by analogy with the parity basis method, we also act with a x on all qubits in 
the appropriate U(j); that is, we apply the operator Xjj(j)- This has the effect of updating 
all the qubits that store a set of orbitals including orbital j. The size of U(j) also scales like 
O(logn) [13]. To summarize: to represent aj or aj in the Bravyi-Kitaev basis, for j even, 
we act with a z on all qubits in P(j), on qubit j, and with a x on all qubits in U(j): 

a] = X u{j) <g) Qf <g> Z P{j) = \{X u{j) ® Xj <g> Z P {j) - iXuQ) <8> Yj <g> Z Pij) ); (32) 
aj = Xu(j) <E> Qj ® Z P Q) = \{Xu(j) ® Xj ® Z P {j) + iXu(j) <8> Yj (g) Z P ^). (33) 
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In the next section, we will consider the case for which j is odd. 

B. Representing a~p in the Bravyi-Kitaev basis for j odd 

To represent the creation or annihilation of a particle in orbital j in the Bravyi-Kitaev 
basis, for j even, we could simply act with Q ± on qubit j because that qubit stores only the 
occupation of orbital j. For j odd, qubit j stores a partial sum of occupation numbers of 
orbitals including, but not limited to, orbital j. Thus, in this case the state of Bravyi-Kitaev 
qubit j is either equal to the occupation of orbital j (if the parity of the other orbitals that 
it stores is even), or opposite to that of orbital j (if the parity of the other orbitals that it 
stores is 1). Thus, whether representing the creation or annihilation of a particle in orbital 
j requires that we act with Q + or Q~ on qubit j in the Bravyi-Kitaev basis depends on the 
parity of all occupation numbers other than fj that are included in the partial sum bj — i.e. 
the parity of the flip set of index j. If the parity of the set of qubits with indices in F(j) is 
even, then the creation or annihilation of a particle in orbital j requires acting with Q + or 
Q~ , respectively, as usual. But if the parity of this set of qubits is odd, then the creation of 
a particle requires acting with Q~ and the annihilation of a particle requires acting with Q + . 
The Bravyi-Kitaev analogues to the qubit creation and annihilation operators are therefore: 

flf = Qf ® E m - QJ ® O m = 1 2 {X 3 ® Z m =F iYj). (34) 

The updating procedure in this case in which j is odd works in exactly the same way as 
it does in the case that j is even. In applying the parity operator, however, we need only 
consider the qubits that are in P(j) but not in F(j), because the relative sign in the flf 
operator implicitly calculates the parity of the subset of the parity set that is also in the flip 
set of index j. It is convenient to therefore introduce the new "remainder set": 

R(j) = P(j)\F(j). (35) 

Thus, the fermionic creation and annihilation operators acting on orbital j for j odd are 
represented in the Bravyi-Kitaev basis as follows: 

a] ee X V (j) <g) 11+ <g) Zrq) = \{X V {j) <g> Xj <g) Z P (j) - iX V {j) <g> Yj ® Z r q)); (36) 
a,j ee X u{j) ® f[~ ® Z R(J ) = \{X u{j) <g> Xj ® Z P (j) + iX u{j) ® Yj ® Z R{j) ). (37) 
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It is evident by inspection that the only difference in the algebraic form of the operators 
between the even- and odd-indexed cases is that the second term involves Zpu\ for the even 
case, but Z R ^ for the odd case. Therefore we define: 

PU)={ P<J) lUiSeVen; (38) 
R(j) if j is odd. 

Now the fermionic creation and annihilation operators acting on arbitrary j are represented 
in the Bravyi-Kitaev basis as: 



a] = Xu(j) ® fl+ <g> Z m = \{X u{j) <g> Xj <g> Z P{j) - iX u(j) <g> Yj <g> Z p{j) ); (39) 
aj = X u{3) ® ftj ® Z m = \{X u{j) ® X, ® + zX y(j -) <8) Yj ® Z p{j) ). (40) 

These are useful basic results, but the operators that appear in the molecular electronic 
Hamiltonian are actually products of these creation and annihilation operators. In the next 
section, we derive general expressions for products of these second-quantized operators. 

VI. PAULI REPRESENTATIONS OF SECOND-QUANTIZED OPERATORS IN 

THE BRAVYI-KITAEV BASIS 

In this Section we derive simplified algebraic expressions for classes of Hermitian second- 
quantized fermionic operators in the Bravyi-Kitaev basis. The five relevant classes of oper- 
ators are summarized in Table |TJ We will give complete compact algebraic expressions for 
only the number operators and the Coulomb and exchange operators. It is not possible 
to give the algebraic form for the remaining three classes of operators without considering 
an impractical number of sub-cases, so we opt to give general expressions for products of 
the form a\aj, and show how to use these results to generate algebraic expressions for the 
remaining classes of operators. 

A. Number operators: ha a|oj 

The number operators are of the form ha a\di and have eigenvalues corresponding to the 
occupation number of orbital i. We would like to find a simplified expression for this class 
of operators in the Bravyi-Kitaev basis. 
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Operator 


Second quantized form 


Number operator 




Coulomb/exchange operators 


j j t j j 


Excitation operator 


(ajcij + ajaj) 


Number-excitation operator 




Double excitation operator 


hijki (ol a i afca/ + a \ a \ a j a i) 



TABLE I: The five classes of Hermitian second quantized operators that appear in electronic 
Hamiltonians. In general the overlap integrals hij and hij^i may be complex. 

Given the results of Section [Vj we can write the following: 

a\ai = \{X u{i) ®Xi® Z P[i) - iX u{i) ®Y l ® Z p[i) ) (41) 
x \{X U({) ®Xi® Z P{i) + iX u(i) ®Yi® Z p(i) ). 

Given that a x a x = o^cyy = a z a z = 1, it follows that (X 5 ) 2 = (Y s ) 2 = {Z s f = 1. We are 
left with: 

a\ ai = i[l + tiXiYi) ® Zp (i) \ p{i) - i(Y t Xi) ® Z P(i) \ p{€} + 1] (42) 

= \{1- Z t ®Z P{i)Xp{i) ). (43) 

Now, when i is even, p{i) = P{i), and so P{i) \ p{i) = 0. When % is odd, p{%) = R(i), and 
so P(i) \ p(i) = F(i). Conveniently, F(i) = for i even, so if we define the following: 

Fit) = F{i) U {i}, (44) 

then we can represent the number operators for arbitrary i (even or odd) as follows: 

a Ui = \( 1 ~ z m)- ( 45 ) 

In the next section we consider the Coulomb and exchange operators. 

B. Coulomb and exchange operators: hijji a\ajajai 

The Coulomb operators are of the form while the exchange operators are of 

the form Since these two kinds of operators can be grouped together 
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algebraically, we consider them as one case. The fermionic anti commutation relations ensure 
that i)(aj(ij). Thus, we can consider the Coulomb and exchange 

operators as a product of two number operators. With the result from Section VIA we can 
write the following: 



ala]ajai = §(1 - Z F{i) ) x |(1 - Z F(J) ) (46) 

= K 1 ~ Z F{i) ~ Z *W + Z F(i) Z F(i))- ( 47 ) 

Any overlap between supp(Zp(j)) and supp(Z F ^)), where supp(O) is the support of the 
operator O, i.e. those tensor factors on which it acts nontrivially, will result in the local 
product a z a z = 1. Thus, we only actually need to act with a z on the union of F(i) and F(j) 
minus their intersection, i.e. the symmetric difference of these two sets. Thus we define the 
following notation: 

^^F(0AF(j) = (i^)U F(j)) \ H F(j)) . (48) 
We can then give the algebraic expression for the Coulomb and exchange operators: 

a\a]ajai = -(1 - Z F(i) - Z F(J) + Zp^). (49) 
In the next section we consider general products of the form a\aj. 



C. Products of the form a]aj 

We can assume without loss of generality that i < j. The algebraic form for products 
of this kind depends on the parity of the indices. There are four cases and we will work 
through the first case in detail, and simply present the results for the other cases. 

Using the result of Section |Vl we obtain the following when i and j are even: 



a] a 



i^j — 2 



\ (X u{i) ®Xi® Z P(i) - iX u(i) ®Yi® Z P{i) ) (50) 



x \{X uf j) <g> Xj <g> Z P( j) + iX u{j) ® Yj <g> Z P (j)). 



For each of the four terms resulting from multiplying out the operators in equation (50) 
above, we must consider what products of local qubit operators can result. There are three 
potential sources of local qubit operator products: overlap between the update set of qubit 
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i and the update set of qubit j, overlap between the update set of qubit i and the parity set 
of qubit j, and overlap between the parity set of qubit i and the parity set of qubit j. Any 
overlap between the update sets of qubits % and j will result in the local product a x a x = 1; 
any overlap between the update set of qubit % and parity set of qubit j will result in the 
local product ±ia y ; and any overlap in the parity sets of qubits i and j will result in the 
local product a z a z = 1. Thus we define the following sets: 

£% = [/(z)A[/(j), « 13 E^)n?(j), P°=P(i)AP{j). (51) 

Note that in the case that i and j are even, we do not need to consider the possibility that 
j G U(i) because U(i) contains only odd elements. Similarly, we do not need to consider the 
possibility that i G P(j), because P(j) for j even contains only odd elements. 

As an example, we will show how to use the sets defined above to simplify the term 
{Xuu) <S> Xi <S> Zp(q)(Xu(j) <S> Xj ® Z P rj\). For this term, we need only apply a x to the set of 
qubits Uij \ ctij U {i,j}, cr y to the qubit with index in (which set in general has at most 
1 element, and in the case that % and j are even always contains 1 element), and a z to the 
qubits in the set PP. \ a, L y Thus, this term simplifies to: 

(Xu(i) ® Xi (g) Z P (i))(Xu(j) <8> Xj ® Z P (j)) = -i Xu^a^ijyYa^Zpo^.. (52) 

Using the same reasoning for the other terms, we arrive at the following result: 

a\a 3 = ±X UiAaij Y aij Zpo M . [Y^ - Xtf - liX.X, + Y d Y t )\. (53) 

This is our result for the case that % and j are even. The algebraic expressions for the other 
cases can be derived in the same manner, with the added complication that the expression 
for the product a\aj varies, depending on if i e P{j) and/or j G U(i). This complication 
results in a proliferation of sub-cases: two for the case that % is odd and j is even, three for 
the case that i is even and j is odd, and four for the case that % and j are odd. The only 
additional sets we need to define are the analogs of P?- for when one or both of the indices 
are odd: 

Pi = P(t) AR(j), P^ = R(i)AP(j), P^ = R(i)AR(j). (54) 

The results for all cases are summarized below in Table [TTJ In the following sub-sections we 
show how to use the contents of Table [TT] to generate algebraic expressions for the excitation 
operators, the number-excitation operators, and the double-excitation operators. 
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D. Excitation operators: hij {a\aj + a^ai) 
Providing for the possibility that the integral hij is complex, we can write: 

hij {a\dj + ajctj) = R{hij}(a\aj + ajaj) + ^s{hij}(a\aj — a]aj). (55) 
Applying this to the case when i and j are even, we find the following: 




+ ^{h^iXjXi + YjYi)}. 

Similar expressions for other cases are easily generated by taking the appropriate form of 
a\dj from Table |nj 

E. Number-excitation operators: hyj k (a\ajaja k + a\a^ajai) 
Due to the fermionic anti-commutation relations, the following is true: 

alajCijCtk + a k aj(ijai = {a\a k + a k ai)(a^aj). (57) 

We see that this is simply a product of an excitation operator and a number operator. We 
have previously given algebraic expressions for both of these classes of operators, so it is 
not difficult to combine them for an expression for the number-excitation operators. Let us 
consider the example when i and k are even. Then we have the following: 

hijjk (di^k + a k a i) a j a j = \X{j ik \a lk Ya ik ZpC [^t{h i jjk}(Yj t X i — Xfcl^) (58) 

+%{hijj k }{X k Xi + Y k Yi)] x |(1 - Z m ). 

To simplify, all we need to consider is the intersection between F(j) and the support of 
(a\a k + a\di). In this case the support of the excitation operator is Ui k U a ik U Pf k U {i, k}. 
The form of the simplification will vary depending on these sets, but the process of reducing 
local operator products by exploiting the relationship between the three Pauli matrices is 
unchanged. In the cases when i and k are not both even, all that changes is the form of the 
excitation operator from Table [IT] that must be used. 
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F. Double-excitation operators: hijki (a\a}jakai + a^a^ajai) 

The double-excitation operators involve four distinct indices, and are obviously the most 
algebraically complicated class of operators we are considering. The impractical number 
of sub-cases depending on the specific combination of indices i,j,k,l means that we only 
outline the procedure for deriving algebraic expressions for this class of operators. The 
fermionic commutation relations ensure that the following is true: 

(alajdkai + a\a\aj(Xi) = (a\ai)(a^ak) + (a\ai)(a} k aj). (59) 

Allowing for the integral hijki to be complex, we can write: 

hijki (a\a)a k ai + a\a\ajai) = [ ^{hijki}(a\aia]a k + aja^a,) (60) 

+ Q{h ijM } {a\aia]a k - aja^a,)]. 

Since (a\aia^a k y = ajaia k aj, we can simply consider the algebraic expression for the product 
of two operators of the form a\aj as given in Table [nj and then add or subtract it to its 
Hermitian conjugate. Each of the operators a\ai and a^a k will fit into one of the ten cases 
presented in Table |TTJ In multiplying out the algebraic expressions for these two products, 
what is important is the set {supp(ajaz) fl supp(a]afc)}. Any qubits in this set will have a 
product of local operators acting on it which must be simplified. 



VII. THE MOLECULAR ELECTRONIC HAMILTONIAN FOR THE 
HYDROGEN MOLECULE IN THE BRAVYI-KITAEV BASIS 

The molecular electronic Hamiltonian ^ may be divided into one and two-electron terms: 
H = hijajaj + - Y h ijkl a\a]a k ai = H {1) + H {2) . (61) 

i,j i,j,k,l 

We treat molecular hydrogen in a minimal basis, so the sums above run over the four 
spin orbitals defined above. These spin orbitals will be indexed through 3, as will be the 
fermionic creation and annihilation operators. We derive the simplified expressions for the 
individual terms of this Hamiltonian in the Bravyi-Kitaev basis. The overlap integrals hij 
and h^ki for < i < 3 are given in Table III These are the same as were used in [11 J and 



were calculated using a restricted Hartree-Fock calculation in the PyQuante quantum chem- 
istry package [32]. With these integrals and the algebraic expressions for second quantized 
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Index parity 




Conditions 




Algebraic expression for a\aj 




i £ P(j) 


j G 17(0 


\ a ij\ 




i,j even 


No 


No 


1 


jXtj .\ a .. Y a .. Z P o\ \YjXi — XjYi — iiXjXi + Yj-y)] 




No 


No 


1 


Y ay [(Y.-JQ - /.Y ; .Y,j Z P o - {X^i + sY^) Z Pf ] 


i odd, j even 












Yes 


No 





\X Vi . Z % [(YjYi - iXjXiYi) Z p0 + (XjXt + iYjX,) Z p2 } 




No 


No 


1 


ix^v^ [-(^-y + iXjXi) z P% + (y,-^ - iYjYi) z P i] 


i even, j odd 


No 


Yes 


1 


\x UiAj [-X Qi .(y - *y.) y Qij z P o w . + (;y - x t ) z^] 




Yes 


Yes 





\X UiAo [(X, - iYi) + (iYi - Xi) Z P i u ,] 




No 


No 


1 


■aXtj. .\fy, . Y OL .-Z(x i A f — iXjXiZpo -\-YjX{Zpi — XjYlZp2 — iYjY{Zpz\ 




Yes 


No 





\x v .. ZiK-iXjYiZpo + y^-y Zpi ) + x j x i z P 2 + iYjXiZps ] 

ij >3 ij 


i,j odd 












No 


Yes 


1 


I^ijV ["^a^ (y -^p2 +iXiZpo)y ay Z aj;/ - {XiZpi - iYiZ P z)Zj\ 




Yes 


Yes 





\Xu Xj [Zii-iYiZpo + XiZ P 2) + Zji-XiZpi +iY l Z pa )] 



TABLE II: The algebraic expressions for general products of the form a]aj in the Bravyi-Kitaev 
basis. These expressions vary in form depending on the parity of the indices i and j, as well as on 
the overlaps between the parity and update sets of the indices. The notation Os is shorthand to 
indicate that the operator O does not operate on the qubits in the set S (i.e. Z P o Zj = Z P o \ ■). 

ii ij w 



operators given in Section VI, we can express the molecular electronic Hamiltonian for H 2 
as a sum of products of Pauli matrices. In the next two subsections we consider the one- 
and two-electron Hamiltonians separately. 
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Integrals 


Value (a.u.) 


^-00 = ^11 


-1.252477 


h 2 2 = h 33 


-0.475934 


^ono = /iiooi 


0.674493 


^2332 = ^3223 


0.697397 


^0220 = ^0330 = ^1221 = ^1331 
= ^2002 = ^3003 = ^2112 = ^3113 


0.663472 


^0202 = ^1313 = ^-2130 = ^2310 = ^0312 = ^0132 


0.181287 



TABLE III: The overlap integrals for molecular hydrogen in a minimal basis. The integrals were ob- 
tained through a restricted Hartree-Fock calculation in the PyQuante quantum chemistry package 
at an internuclear separation of 1.401000 atomic units (7.414 x 10~ n m). 

A. The Bravyi-Kitaev Pauli representation of 

We can write the one-electron terms in the Hamiltonian as: 

= h 00 ala + hua\ai + /i 2 2a-2 a 2 + h 33 a\a 3 . (62) 

Using the expressions for number operators derived in Section |Vj we know that in the 
Bravyi-Kitaev basis, these operators are: 



(63) 
(64) 
(65) 
(66) 



a\a = |(1 - <rg); 
a\ax = |(1 - crfa$); 

a\a 2 = |(1 - af); 
a\a 3 = \{\-alo z 2 a\). 

We now proceed to the simulation of H^ 2 \ 

B. The Bravyi-Kitaev Pauli representation of 

Following the work of Whitfield et al. [llj, H^ 2 ' simplifies to the following expression for 
molecular hydrogen in a minimal basis: 



# (2) = Vl0<44 a l a + ^23324 a 3 a 3 a 2 +^033o44 a 3«0 + /il22ial a 2 a 2CH (67) 
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+ (^-0220 — ^0202)44 a 2«0 + (^1331 ~~ ^1313) 44 a 3 a l + ^-0132 ( a a l a 3 a 2 + 44 a l a o) 

+h^i 2 (a\a\a 1 a 2 +4a|a 3 a ). 
This term in the Hamiltonian is made up of six Coulomb/exchange operators and two 



double-excitation operators. Using Section VI, it is easy to give algebraic expressions for 
the Coulomb and exchange operators: 

44^0 = \(1 - a z Q - afa z + of); (68) 
a\ala 3 a 2 = f (1 - erf - erf erf a\ + a z cx(); (69) 

44a 3 a = |(1 - Crg - cr^C^ + cr|cr|(Tfcrg); (70) 

44a 2 ai = |(1 - erf - a z a^ + (T z afa^); (71) 

44^200 = f(l -erf -crg + erferg); (72) 

44a 3 ai = |(1 - o-fcrfo-i - o-?o"o + ■ ( 73 ) 

The two double-excitation operators are somewhat more complicated. As an example, we 
will derive the Pauli representation of /io3i2(44 a i a 2 + 44 a 3 a o)- Following in Section 
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we consider 44 a i fl 2 as (4 a 2)(4 a i)) a product of two operators of the form a\aj. The term 
4<^2 is of the type when % and j are both even, while the term 4°3 is °f the type when i 
and j are odd, and i G P(j), j G U(i), and \a^\ = 0. Using the appropriate expressions 
from Table [TTJ we find the following: 

4°2 = \((T2<Tl<7o ~ ^l^O - wZOlOo ~ ^2^0 ) 5 ( 74 ) 

4 o 3 = K-icrfo'lO'o + o&l - ^l^o + ^M)- (75) 

Now we note that supp(4&2) H supp(a|a 3 ) = {2, 1,0}, and so we must expect to simplify 
local operator products on qubits with these indices. Taking the product, we find the 
following: 

ala 2 a{a 3 = — ( afcrg - zof 0$ + crfo-fc^ - ia^afa^ (76) 
16 

+ ia\a% + af cr^ + icr%a z cr$ + erfcrfc^ 

+ Cftfftfo - i°iO2°0 + °'3 Cr 2 Cr l O '0 - 

+ i(T%(T%(T% + ofer]^ + wfcrfcrfcrg' + erf <jf ofc^) . 

Since the integral /10132 is real, we can simply add the above result to its Hermitian 
conjugate to find the expression for the double-excitation operator. Repeating the above 
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procedure for the second double excitation operator, we arrive at the following results: 

44°l a 2 + 44^0 = o( -O-f°"0 + cr 2°"l cr _ a 2°l + ^2^0 ~ °l°2 a l ( 77 ) 

o 

44«302 + 44 a l a = g( 0%Oq + ^2 a l a + cr 2 cr + ^2^0 +^2^0 ( 78 ) 



Thus, using the integrals from Table III and the Pauli expressions for the number opera- 
tors derived in Section VII A, as well as the Coulomb/exchange operators and the double- 
excitation operators derived in this section, we can represent the molecular electronic Hamil- 
tonian for the hydrogen molecule as a sum of products of Pauli matrices in the Bravyi-Kitaev 
basis: 

H BK = -0.81261 1 + 0.171201 a z + 0.16862325 o\ - 0.2227965 o* + 0.171201 a\a% 

+0.12054625 erf <rg + 0.17434925 a\a\ + 0.04532175 a%a z a% + 0.04532175 o%afo% 
+0.165868 o-faf erg + 0.12054625 0-fo-fcr* - 0.2227965 affiaf 
+0.04532175 (rfof of (j* + 0.04532175 o^ofag + 0.165868 a'a'affi. (79) 

This Hamiltonian is isospectral to the Jordan- Wigner derived Hamiltonian [TT] : 

H JW = -0.81261 1 + 0.171201 a z + 0.171201 a\ - 0.2227965 a z - 0.2227965 o% 
+0.16862325 afa z + 0.12054625 a z 2 a z + 0.165868 erf a\ + 0.165868 a z 3 a z 

+0.12054625 crfcrf + 0.17434925 <j%o% - 0.04532175 o$o%ofa% 
+0.04532175 of ala\o% + 0.04532175 a y 3 a^a y - 0.04532175 of afcf <7*- (80) 



Writing the electronic Hamiltonians in the form of equations (79) and (80) allows for a com- 
parison of the computational resources required to simulate them on a quantum computer. 
Not all tensor products of Pauli matrices that appear in these Hamiltonians commute with 
one another, so exponentiating them requires the use of a Trotter approximation. The next 
section details the Trotterization process for the Hamiltonian in the Bravyi-Kitaev basis. 

VIII. TROTTERIZATION 

Ideally, one could simulate the propagator e~ lHt , where H = J2k hk, by sequentially 
exponentiating the individual terms hk on a quantum simulator. However, e~ lHt = T7 e~ thkt 
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only in the case that the set of hk all mutually commute. Both the Bravyi-Kitaev and 
Jordan- Wigner Hamiltonians contain terms that do not commute with one another, and 
so a Suzuki- Trotter approximation must be used. The first four orders of Suzuki- Trotter 
formulae are |27j : 

e (A+B)t _ ^At/n e Bt/ny + Q^ty ( gl ) 
e (A+B)t ^ ( e At/2n e Bt/n e At/2ny + (t(At) 2 ); (82) 
e (A+B)t ^ ( e ^At/n e lBt/n e lAt/n e ^Bt/n e ^At/n e Bt/ny + (t(Atf)- (83) 
e (A+B)i ^ (J| ^gwAt/angRBt/ngMt/^n + 0(t(At) 4 ), (84) 
where in the 4th order equation, the constants are given by: 

Pi = P2 = Pi = P5 = 4 _ 1 4l/3 , p 3 = l-4pi. (85) 

The terms of both the Bravyi-Kitaev Hamiltonian and the Jordan- Wigner Hamiltonian can 
be broken into two subsets, where the terms in each subset all mutually commute but the 
subsets do not commute with one another. These groups are as follows: 

H BK ,z= - 0.81261 1 + 0.171201 + 0.16862325 a\ - 0.2227965 a z + 0.171201 a\a% 
+ 0.12054625 a z 2 a z + 0.17434925 afo-f + 0.165868 o z 2 a\al 

+ 0.12054625 o-|o-|o-g - 0.2227965 o-|o-|o-f + 0.165868 o-fo-fo-fo-g; (86) 

H B k,xy = 0.04532175 ofofofi + 0.04532175 o\o\ol + 0.04532175 affiafof; 

+ 0.04532175 crfo-fo^; (87) 

H JWZ = - 0.81261 1 + 0.171201 a z G + 0.171201 o\ - 0.2227965 a% - 0.2227965 erf 
+ 0.16862325 afa z +0.12054625(7^ + 0.165868 a z 2 a z x + 0.165868 
+ 0.12054625 erf of +0.17434925 a z 3 a z ; (88) 



Hj WiXY = - 0.04532175 exf a^ofo^ + 0.04532175 of o-fof of + 0.04532175 a%o%ofo% 

0.04532175 of crfof of . (89) 

To understand what computational resources are required for exponentiating operators of 
this kind, consider the example of the exponentiation of a fourfold product of a z matrices, 
e iO z ®cr z (g)CT z ®cr z )^ -which j s depicted in a circuit diagram in Figure [3] [28] . 
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FIG. 3: A demonstration of how to exponentiate tensor products of Pauli matrices. First, the 
parity of the four qubits is computed with CNOT gates, and then a single-qubit phase rotation R z 
is applied. Then, we uncompute the parity with three further CNOT gates. 

In general, an n-fold tensor product of Pauli-Z matrices will require 2(n — 1) CNOT 
gates and one single-qubit gate (SQG) to exponentiate on a quantum computer. If there are 
Pauli-X or -Y matrices in the tensor product, we must apply the single-qubit Hadamard or 
R x gate to change basis to the X or Y basis, respectively, before we compute the parity of 
the set of qubits with CNOT's, and also apply the inverse gates as part of the uncomputing 
stage [2H]- These gates are given by: 



H 



1 

V2 



1 1 
1 -1 



R, 



1 

V2 



1 i 
i 1 



(90) 



Thus, each non-a 2 term in a tensor product of Pauli matrices adds 2 single-qubit gates to 
the cost of exponentiation. For example, the circuit for exponentiating the term erf erf erf erg 
is depicted in Figure [4] 

Using the resource counting methods detailed above, we can count the number of single- 
qubit gates (SQG's) and CNOT gates required to exponentiate (for arbitrary propagation 
time) the subsets of the Hamiltonians for both encodings. The results of this analysis are 
in Table [T3 

We now have the tools to compare the number of gates required to compute the ground 
state eigenvalue of either the Bravyi-Kitaev Hamiltonian or the Jordan- Wigner Hamiltonian 
to chemical precision (±10 -4 a.u). Due to the small size of our model of the hydrogen system, 
it is easy for a classical computer to simulate the behavior of the quantum simulator. The 
true propagator U = e~ lHt can be computed to sufficient precision by a matrix exponential 
function in Mathematica or a similar software package. Time evolution of the ground state 
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FIG. 4: A demonstration of how to exponentiate tensor products of Pauli-X and -Y matrices. 
First, the qubits are put in the correct basis by the application of R x or Hadamard gates. Then, 
the parity of the four qubits is computed with CNOT gates, and then a single-qubit phase rotation 
R z is applied. Then, we uncompute the parity with more CNOT gates, and finally change back to 
the computational (Z) basis. 





SQG's 


CNOT's 


Totals 


H B k,z 


10 


24 


34 


Hbk,xy 


20 


20 


40 


Totals 


30 


44 


74 


Hjw,z 


10 


12 


22 


Hjw,xy 


36 


24 


60 


Totals 


46 


36 


82 



TABLE IV: The number of single-qubit gates and CNOT gates required to exponentiate subsets 
of the electronic Hamiltonian for the hydrogen molecule, represented in terms of spin variables 
through either the Bravyi-Kitaev transformation or the Jordan- Wigner transformation. 



by the true propagator will result in phase evolution: 

U\1> g ) = e-^U>. (91) 
We can therefore compute the exact eigenvalue as follows: 

(^IW ff > = (^ 9 \e- iE ^ g ) = e-^K (92) 

We set the propagation time to unity, and extract the true eigenvalue E g from the complex 
phase e~ %E 3 . To approximate the eigenvalue, we use a Suzuki- Trotter approximation to the 
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true propagator, U, and perform an analogous procedure: 

e"^ 4 . (93) 



\W\^ 9 )\ 

The approximation to the true ground state eigenvalue, E g , becomes better as we increase 
the number of Trotter steps n. Figure [5] below plots the estimated eigenvalues of the minimal 
basis Jordan- Wigner and Bravyi-Kitaev Hamiltonians as a function of the number of gates 
required, for the first four orders of Suzuki- Trotter formulae. 

We now compare this result to previous estimates. The benchmark is the gate count 
given in [11] for approximating the Jordan- Wigner Hamiltonian's ground state eigenvalue. 
It is clear from Figure [5] that our first order approximation requires ~ 900 gates to obtain 
chemical precision for the Jordan- Wigner Hamiltonian, while the gate estimate in [TTJ was 
about 500 for the same task. This discrepancy arises from the fact that any number of 
variants on the first order Suzuki- Trotter formula could have been used in [UJ. Given a 
noncommuting set of Hamiltonian terms, there is some optimal ordering that will produce 
the best accuracy. It is not possible to know in advance which ordering is optimal, and 
given that the number of terms in an electronic Hamiltonian scales as 0(n 4 ), in general it 
is difficult to optimize over the space of possible orderings. We have used the most naive 
variant of the first order Suzuki- Trotter formula in Figure [5] 

e -iHt _ e -i(H z+ H XY )t _ , _i# 2 £ e _»jy xv Av„ 

However, due to the small size of our model of the hydrogen molecule, it is easy to find 
an ordering that produces better accuracy. A second, more sophisticated, variant of the 
first order formula is to arrange the terms in Hz and Hxy i n order of descending coefficient 
magnitude. For example, for the Bravyi-Kitaev Hamiltonian, we have: 

H z : {h Z0 , hzi, h Z 2, ■ ■ ■} = {-0.81261 1, -0.2227965 erf, -0.2227965 ofofof, . . •}; (95) 
Hxy ■■ {hxYo, h XY i, h XY2 , . . .} = {0.04532175 ofo^og, 0.04532175 of o*o£, . . .}. (96) 

Then, we approximate the propagator by alternately exponentiating one term from the 
ordered list of Hz terms and one term from the ordered list of Hxy terms until we have 
used all terms from Hxy- Then we exponentiate the rest of Hz' 

e -iHt ^ f e -ihzo± e -ih X Yo± e -ihzi± e -ih X Yi± . . . e -ih X Y3^ e -ih Z i^ e -ih zs ± . . 
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FIG. 5: The approximation to the ground state eigenvalue, for both the Bravyi-Kitaev Hamiltonian 
(squares) and Jordan- Wigner Hamiltonian (circles), as a function of the number of gates required. 
The solid curves are the first order Suzuki- Trotter approximations, the dot-dashed second order, 
the dotted third order, and the dashed fourth. The dotted horizontal line represents the true 
eigenvalue, while the solid lines above and below represent the bounds for chemical precision. 

With this method, we find that the number of gates required to obtain a chemical precision 
estimate of the ground state eigenvalue of the Jordan- Wigner Hamiltonian is 300, fewer 
than the result from [TT] . Figure [6] compares the eigenvalue approximations for the naive 
first order method and the more sophisticated variant. 

The point is that the systematic advantage of the Bravyi-Kitaev method over the Jordan- 
Wigner method is not obscured by the kind of term-ordering optimization that we have 
demonstrated above. Exponentiating the Bravyi-Kitaev Hamiltonian requires 74 gates per 
first order Trotter step (of any variant), while the Jordan- Wigner Hamiltonian requires 82 
gates per first order Trotter step. To obtain a precision of ±10 -4 a.u to the true eigenvalue 
with the naive first order Suzuki- Trotter approximation requires 11 Trotter steps for both 
the Bravyi-Kitaev and Jordan- Wigner Hamiltonian, for a total cost of 814 gates versus 
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FIG. 6: The approximation to the ground state eigenvalue, for both the Bravyi-Kitaev Hamiltonian 
(squares) and Jordan- Wigner Hamiltonian (circles), as a function of the number of gates required. 
The solid curve is the naive first order Suzuki- Trotter approximation, while the dashed curve is 
the result from alternating the noncommuting terms. The dotted horizontal line represents the 
true eigenvalue, while the solid lines above and below represent the bounds for chemical precision. 
The ground state eigenvalue of the Bravyi-Kitaev Hamiltonian can be approximated to chemical 
precision with 222 gates, while it takes 328 gates to do the same for the Jordan- Wigner Hamiltonian. 



902 gates. With the noncommuting terms intermixed, it takes only 3 Trotter steps to 
obtain the same precision for the Bravyi-Kitaev Hamiltonian, and 4 Trotter steps for the 
Jordan- Wigner Hamiltonian. Thus, if we intermix the noncommuting terms, the Bravyi- 
Kitaev transformation allows one to utilize 222 gates instead of the 328 gates required by 
the Jordan- Wigner transformation to obtain an equally precise estimate of the hydrogen 
molecule's ground state eigenvalue when using a first order Suzuki- Trotter approximation. 
When using higher-order Suzuki- Trotter approximations to obtain better than chemical 
precision, the gate savings increases (Fig. W\. 
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FIG. 7: The gate savings of using the Bravyi-Kitaev method instead of the Jordan- Wigner method, 
as a function of the precision in the estimate of the ground state eigenvalue for the first four orders 
of Suzuki- Trotter formulae. The vertical line is the threshold error for chemical precision. The 
triangle data points are first order, the squares second, the circles third, and the diamonds fourth. 

IX. CONCLUSIONS 

In this paper we have worked out a detailed application of the Bravyi-Kitaev transforma- 
tion to Hermitian second quantized operators that appear in quantum chemical Hamiltoni- 
ans. We suggest that this transformation should replace the Jordan- Wigner transformation 
for fermionic quantum simulation algorithms. We have demonstrated that the Bravyi-Kitaev 
transformation results in a small reduction in the number of gates, from 328 gates to 222 
gates, required to implement a quantum simulation algorithm for electron dynamics in the 
simplest possible molecular system of H2 in a minimal basis. 

In some sense, molecular hydrogen in a minimal basis is a poor showcase of the power of 
the Bravyi-Kitaev transformation. Our description of this molecule utilizes four molecular 
orbitals, and hence four qubits. The spin Hamiltonians we derive using either the Bravyi- 
Kitaev transformation or the Jordan-Wigner Hamiltonian involve four-local Pauli tensor 
products, the result being that the cost of simulating time evolution under the Bravyi- 
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Kitaev Hamiltonian on a quantum computer is only slightly reduced from that for the 
Jordan- Wigner Hamiltonian. However, were we to use a more sophisticated description of 
the H 2 — for example, with eight molecular orbitals — the Jordan- Wigner spin Hamiltonian 
would contain up to eight-local Pauli tensor products, while the Bravyi-Kitaev spin Hamil- 
tonian would not. Given the asymptotically better O(logn) scaling of the Bravyi-Kitaev 
method as compared to the 0{n) scaling of the Jordan- Wigner transformation, the differ- 
ence between the two methods will become greater for larger basis sets and larger molecules 
- the simulation of which is, after all, is the true goal of quantum simulation for quantum 
chemistry, since the small molecules are within the reach of conventional computers. How- 
ever, by showing that the Bravyi-Kitaev method is more efficient for the smallest conceivable 
chemical system, we have demonstrated that there is no algorithmic overhead inherent to the 
Bravyi-Kitaev method that must be overcome by scaling up the size of problems to which 
it is applied. We have demonstrated the superior efficiency of the Bravyi-Kitaev trans- 
formation for all quantum chemical simulations. Thus, making use of the Bravyi-Kitaev 
transformation for fermionic quantum simulation will make simulations of larger molecules 
and with larger basis sets more readily accessible to experiment. 
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